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ABSTRACT 

We study the evolution of the scaling relations between maximum circular velocity, 
stellar mass and optical half-light radius of star-forming disk-dominated galaxies in 
the context of ACDM-based galaxy formation models. Using data from the literature 
combined with new data from the DEEP2 and AEGIS surveys we show that there is a 
consistent picture for the evolution of these scaling relations from z ~ 2 to z — 0, both 
observationally and theoretically. The evolution of the observed stellar scaling rela- 
tions is weaker than that of the virial scaling relations of dark matter haloes, which can 
be reproduced, both qualitatively and quantitatively, with a simple, cosmologically- 
motivated model for disk evolution inside growing Navarro-Frenk- White dark matter 
haloes. In this model optical half-light radii are smaller, both at fixed stellar mass 
and maximum circular velocity, at higher redshifts. This model also predicts that the 
scaling relations between baryonic quantities (baryonic mass, baryonic half-mass radii, 
and maximum circular velocity) evolve even more weakly than the corresponding stel- 
lar relations. We emphasize, though, that this weak evolution does not imply that 
individual galaxies evolve weakly. On the contrary, individual galaxies grow strongly 
in mass, size and velocity, but in such a way that they move largely along the scaling 
relations. Finally, recent observations have claimed surprisingly large sizes for a num- 
ber of star-forming disk galaxies at z ~ 2, which has caused some authors to suggest 
that high redshift disk galaxies have abnormally high spin parameters. However, we 
argue that the disk scale lengths in question have been systematically overestimated 
by a factor ~ 2, and that there is an offset of a factor ~ 1.4 between Ha sizes and 
optical sizes. Taking these effects into account, there is no indication that star forming 
galaxies at high redshifts (z ~ 2) have abnormally high spin parameters. 

Key words: galaxies: evolution - galaxies: formation - galaxies: fundamental pa- 
rameters - galaxies: haloes - galaxies: high redshift - galaxies: spiral 



1 INTRODUCTION 



* dutton@uvic.ca 

t CITA National Fellow 

$ Hubble Fellow 

§ Spitzcr Fellow 



In the current paradigm of galaxy formation, galaxy disks 
are considered to form from the accretion of gas inside hier- 
archically growing cold dark matter (CDM) haloes (White & 
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Rees 1978). The dark matter (DM) and gas acquire angular 
momentum via tidal torques in the early universe (Peebles 
1969). When the gas accretes onto the central galaxy, this 
angular momentum may eventually halt the collapse and 
lead to the formation of a rotationally supported disk (Fall 
& Efstathiou 1980). Under the assumption that the total 
specific angular momentum of the precollapse gas is simi- 
lar to that of the DM and is conserved during collapse, this 
picture leads to predictions of present-day disk sizes that 
are in reasonable agreement with observations (Blumenthal 
et al. 1984; Dalcanton et al. 1997; Mo, Mao, & White 1998; 
Firmani & Avila-Reese 2000; de Jong & Lacey 2000; Pizagno 
et al. 2005; Dutton et al. 2007). 

The observational relation between the sizes of galaxy 
disks and their characteristic rotation velocities (the RV 
relation) provides a measure of the specific angular momen- 
tum of galaxy disks (e.g., Navarro & Steinmetz 2000). In the 
simplest model for disk galaxy evolution, the sizes and rota- 
tion velocities of galaxy disks are proportional to the sizes 
and circular velocities of their host dark matter haloes, and 
the evolution of disk sizes at fixed circular velocity scales in- 
versely with the Hubble parameter (e.g., Mo, Mao, & White 
1998). For a ACDM cosmology (with Q M = 0.3, Q A = 0.7) 
this simple model predicts that disks should be a factor of 
~ 1.8 smaller, at fixed circular velocity, at z = 1, compared 
to z = 0, and a factor of ~ 3.0 smaller at z = 2. For a stan- 
dard CDM cosmology (with Qm = 1, = 0) the predicted 
evolution is even stronger: a factor of ~ 2.8 to z = 1 and a 
factor of ~ 5.2 to * = 2. 

Early observations indicated evolution in the RV rela- 
tion with a scaling of ~ (l + z) -1 out to z ~ 1, in agreement 
with theoretical expectations (Mao, Mo & White 1998). 
More recently, Bouche et al. (2007) found that the disk 
scale size-rotation velocity relation at z ~ 2.2 has roughly 
the same zero point as at z — 0, in disagreement with sim- 
ple theoretical expectations. While there are concerns about 
sample selection, small statistics, errors at high-z, etc., it is 
nonetheless appropriate to start thinking about the implica- 
tions of this result for galaxy formation. A proposed expla- 
nation for this non-evolution is that disk galaxies at z ~ 2.2 
have spin parameter a factor of 3 higher than those pre- 
dicted by ACDM (Bouche et al. 2007; Burkert et al. 2009). If 
confirmed, this could imply a significant modification of the 
current paradigm of galaxy disk formation, either through a 
decoupling of the angular momentum of baryons and dark 
matter, or more drastically by changing the mechanism by 
which galactic angular momentum is generated. 

A decoupling between the angular momentum of 
baryons and dark matter haloes is actually seen in cos- 
mological hydrodynamical simulations of disk galaxy for- 
mation (e.g., Navarro & Steinmetz 2000; Sales et al. 2009; 
Piontek & Steinmetz 2009). However, the result is that 
baryons lose angular momentum to the halo, which only 
makes high disk spin parameters harder to explain. A pos- 
sible mechanism for increasing the specific angular momen- 
tum of disk galaxies is the removal of material with low 
specific angular momentum though galactic outflows. This 
mechanism has been shown theoretically to be effective in 



1 The spin parameter is a dimensionless measure of the specific 
angular momentum of dark matter halos. 



dwarf galaxies (Mailer & Dekel 2002; Dutton 2009; Gover- 
nato et al. 2010). Galactic outflows appear ubiquitous in star 
forming galaxies at redshift z ~ 3 (Shapley et al. 2003) and 
2 = 1.4 (Weiner et al. 2009). But whether they are powerful 
enough to increase the specific angular momentum of mas- 
sive (Vrot = 200 kms -1 ) galaxy disks by a whole factor of 3 
remains to be seen. 

The relation between galaxy size and stellar mass (the 
RM relation) of disk-dominated galaxies also shows only 
weak evolution to z = 1 (Barden et al. 2005). This weak 
evolution can be understood in the ACDM scenario as a 
consequence of the increase of dark halo concentrations since 
z — 1 (Somerville et al. 2008) and also because individual 
disk galaxies tend to evolve roughly along the RM relation 
(Firmani & Avila-Reese 2009). At higher redshifts, there is 
evidence of a factor of ~ 2 decrease in half-light sizes of 
disk-dominated galaxies at fixed stellar mass between z = 
and z ~ 2.5 (Trujillo et al. 2006). Again, this evolution is 
weaker than the evolution in the halo RM relation, but it 
is in agreement with the models of Somerville et al. (2008) 
and Firmani & Avila-Reese (2009). We note that elliptical 
galaxies have been found to evolve even more strongly with 
redshift (e.g., Trujillo et al. 2006; van Dokkum et al. 2008), 
and thus there is clear evidence that galaxies of all types 
were smaller, at fixed stellar mass, at higher redshifts. 

It is possible that the non-evolution of the RV relation 
and evolution of the RM relation of disk-galaxies are consis- 
tent. This would require evolution in the MV relation (i.e., 
the stellar mass Tully- Fisher relation, Tully & Fisher 1977), 
with higher M, at fixed V, at higher z. However, the oppo- 
site has been claimed, with ~ 0.4 dex lower M sta r at fixed 
Vmax at z ~ 2.2 compared to z — (Cresci et al. 2009). 
Thus, the various data sets for the evolution of the VMR 
relations are inconsistent. 

In this paper, we examine the VMR relations of disk 
galaxies using data from the literature as well as new results 
from the Deep Extragalactic Evolutionary Probe 2 survey 
(DEEP2, Davis et al. 2003; Newman et al. in prep). We 
discuss, and resolve, sources of discrepancies among the var- 
ious data sets, and compare the data to predictions of a 
simple ACDM-based disk-galaxy evolution model. We show 
that the observations can be reproduced by a model with 
constant spin parameter in agreement with predictions from 
ACDM, and a constant galaxy mass fraction. Hence there 
is no need to invoke abnormally high spin parameters in 
order to explain the scaling relations of disk galaxies at 
2 ~ 2. Throughout we assume a flat ACDM cosmology with 
(n M ,n A ,h) = (0.3,0.7,0.7). 



2 OBSERVATIONS 

2.1 Evolution of the stellar mass Tully-Fisher 
relation 

The relation between stellar mas^] (or luminosity) and rota- 
tion velocity (or linewidth) is also commonly known as the 
Tully-Fisher relation (TF; Tully & Fisher 1977). We define 
the velocity as the maximum rotation velocity, Vmax. Using 

2 Note that stellar mass throughout is total stellar mass, not disk 
stellar mass. 
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Table 1. Evolution of the stellar mass- velocity relation of star 
forming galaxies relative to z = 0.0, using DEEP2 data from 
Kassin et al. (2007). 



rcdshift range 


median z 


Alog 10 M sta r|Sb.5 


N 


0.2-0.5 


0.35 


-0.06 ±0.09 


19 


0.5-0.8 


0.73 


-0.03 ± 0.11 


25 


0.8-1.1 


0.95 


-0.15 ± 0.09 


29 



slit spectroscopy from the DEEP1 survey (Vogt et al. 2005; 
Weiner et al. 2005) , together with optical and A"-band imag- 
ing, Conselice et al. (2005) found that, at fixed Vmax, stellar 
masses are lower by 0.07±0.12 dex at z ~ 0.45, and lower by 
0.11 ±0.13 dex at z ~ 0.85, compared to the z = relation 
from Bell & de Jong (2001). This weak evolution from z ~ 1 
to z ~ was confirmed by Kassin et al. (2007) , who studied 
the stellar mass TF relation for a larger sample of galaxies 
from the DEEP2 survey, but adopting a different velocity 
indicator, 5*0.5 = (0.5VJ.o t ± (J 2 ) 1//2 , which combines ordered 
motions (i.e., rotation, Kot) and disordered motions (i.e., 
dispersion, a) (Weiner et al. 2006a). This new parameter 
allowed inclusion of mergers and incompletely settled galax- 
ies, which is useful at high redshifts were galaxies maybe 
dynamically "young" . 

Here we use data from Kassin et al. (2007) and calcu- 
late the evolution of stellar mass at fixed So. 5 for galaxies 
with So. 5 > 90kms _1 (i.e., log 10 Vmax/Jkms" 1 ] > 2.1) and 
inclinations between 45 and 70 degrees. As a comparison re- 
lation at redshift z = we use the TF relation from Bell 
& de Jong (2001), corrected to a Chabrier (2003) IMF as 
given below in H2.3I and with V to t converted into So. 5 as- 
suming a = 0. We note that since the galaxies at z — 
are rotation dominated, adopting a more realistic value of 
a ~ lOkrns -1 (for the cold atomic hydrogen disk) for these 
galaxies will not change the z — relation by any significant 
amount. The median offsets in stellar mass of the DEEP2 
data with respect to the z = relation, in three redshift 
bins from 2 = 0.2 to z = 1.1, are given in Table [1] These 
results are consistent with no evolution or at most a weak 
decrease in stellar masses at fixed circular velocity at higher 
redshifts. 

A relatively strong evolution in the stellar mass TF rela- 
tion since z ~ 0.6 has been reported by Puech et al. (2010). 
These authors find an increase of 0.34 dex in stellar mass at 
fixed velocity since z ~ 0.6. However, the same authors find 
no evidence for evolution in baryonic mass at fixed veloc- 
ity (the baryonic Tully-Fisher relation) over the same red- 
shift range. The large difference between the evolution of the 
stellar and baryonic TF relations implies that gas fractions 
evolve significantly. However, the inferred gas fractions of 
their galaxies at z ~ 0.6 are on average 30%, which is not 
much higher than that of local galaxies of the same mass. 
This implies that the evolution of the stellar mass TF re- 
lation should be only of order 0.1 dex different than that 
of the baryonic TF relation. The difference in evolution of 
0.4 dex found by Puech et al. (2010) can thus be traced to 
the use of inconsistent local baryonic and stellar mass TF 
relations. 

At even higher redshifts evidence has been reported of 



significant evolution compared to redshift z — 0. Using data 
from SINS (Spectroscopic Imaging survey in the NIR with 
SINFONI, Forster Schreiber et al. 2009), Cresci et al. (2009) 
found that at fixed Vinax, stellar masses at redshift z ~ 2.2 
are lower by 0.41 ± 0.11 dex compared to z = 0. Taken 
together, the observations indicate that there is a modest 
evolution in the zero point of the TF relation out to z ~ 1 
and a stronger evolution from z ~ 1 to z ~ 2. 

2.2 Evolution in the size-stellar mass relation 

Barden et al. (2005) found that there was little or no evo- 
lution in the circularized^] optical half-light size-stellar mass 
relation of disk dominated galaxies (defined as having Sersic 
index n < 2.5) from redshift z ~ 1 to z ~ 0.1. Trujillo 
et al. (2006) measured the evolution of the circularized rest- 
frame V-band half-light radius-stellar mass relation since 
z ~ 2.5, finding that disk-galaxies (Sersic index n < 2.5) 
are a factor of ~ 2 smaller, at fixed M sta r, at z ~ 2.5 than 
at z = 0.1. Williams et al. (2010) measure the evolution of 
the circularized rest-frame /-band half-light radii of galaxies 
with stellar masses greater than 6.3 x 1O 1O M0 from 2 = 0.5 
to z = 2, finding strong evolution for both star-forming and 
non star-forming galaxies. In order to compare to 2 = 0.1 we 
use a mean /-band half-light size of 5.0 kpc. This has been 
determined using our SDSS measurements (see i]2.2.3[l . and 
applying a correction of -0.06 dex to go to rest-frame /-band 
(see 333- 

Here we present new results for the evolution of the 
disk size-stellar mass relation for blue-cloud disk dominated 
galaxies from redshifts 2 = 1. 2 to 2 = 0. 1 using high red- 
shift data from DEEP2 (Davis et al. 2003; Newman et al. , in 
prep) and the All-wavelength Extended Groth Strip Inter- 
national Survey (AEGIS, Davis et al. 2007), and a low red- 
shift comparison sample from the Sloan Digital Sky Survey 
(SDSS, York et al. 2000). The main difference of our study 
with respect to previous studies (e.g. Barden et al. 2005; 
Trujillo et al. 2006) is our use of disk sizes, rather than total 
sizes. We use disk sizes because our main interest in this 
paper is the evolution of galaxy disks. Total sizes, even for 
disk-dominated galaxies, depend on the bulge fraction and 
bulge size, and thus give a measurement that is more diffi- 
cult to interpret. 

2.2.1 DEEP2 data 

Disk sizes have been measured using 2D bulge±disk fits us- 
ing GIm2d (Simard et al. 2002). The bulge component was 
assumed to have a deVaucoulers profile (Sersic n = 4), while 
the disk was assumed to be exponential (i.e., Sersic n = 1). 
The fits were performed simultaneously on F606W (V) and 
F814W (/) single-orbit Hubble Space Telescope (HST) Ad- 
vanced Camera for Surveys (ACS) images. The bulge and 
disk sizes are constrained to be the same in each filter, but 
the bulge and disk fluxes are free to vary. This choice has 
been made to maximize signal to noise, at the expense of 

3 For sizes measured in elliptical apertures, the quoted size, is 
conventionally the major axis size, R. The circularized size is 
then often defined as yJb/aR, where b/a is the minor to major 
axis ratio. 
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using bluer rest-frame wavelengths at higher redshifts: The 
central rest-frame wavelength of the images (i.e., the aver- 
age of the F606W and F814W images) varies from ~ 500 
nm at z — 0.4 to ~ 320 nm at z = 1.2. This may cause 
a small (~ 10%) systematic overestimation of the sizes at 
higher redshifts, because disks tend to have color gradients 
(see q3.3[) . The sizes we present in this paper are major-axis 
(i.e., elliptical aperture) disk half-light sizes, i.e., 1.678 times 
the disk exponential scale length. 

We use both spectroscopic and photometric redshifts 
between 0.2 < z < 1.4. We use high quality ("z-quality" 
3) spectroscopic redshifts from DEEP2 (Davis et al. 2003; 
Newman et al. , in prep), and photometric redshifts based 
on optical to Spitzer/IRAC photometry (J. Huang et al. , 
in preparation). For the spectro-z sample, stellar masses 
were obtained from optical to NIR spectral energy dis- 
tribution (SED) fits from Bundy et al. (2006) which as- 
sumed a Chabrier (2003) initial mass function (IMF). For 
the photo-z sample we use stellar masses calculated using 
a relation between B-band stellar mass-to-light ratio and 
rest-frame (U — B) color, which we calibrate against the 
Bundy et al. (2006) masses. Our calibration is similar to the 
ones reported in Lin et al. (2007) and Weiner et al. (2009), 
but make use of only rest-frame (U — B) colors, rather than 
(U — B) and (B — V). We compute rest-frame (U — B) colors 
and B-band luminosities by applying kcorrect v4.1.4 (Blan- 
ton & Roweis 2007) to DEEP2 B,R,I magnitudes. 

We select blue-cloud galaxies for this study by using the 
"green valley" in the (U — B) color-stellar mass plane as the 
division: 

(U — B) < 1.05 - 0.1(log 10 Af star - 10). (1) 

In addition we select only disk-dominated galaxies with disk 
light fractions greater than 50%. This additional cut only 
removes a small fraction of galaxies because the vast major- 
ity of blue galaxies already have disk fractions greater than 
50%. In order to minimize the effects of extinction on disk 
sizes, colors and stellar masses, while at the same time keep- 
ing a significant sample size, we limit our sample to galaxies 
with disk minor-to-major axis ratios greater than 0.5 (corre- 
sponding to disk inclinations less than 60 degrees, for a zero 
thickness disk). 



2.2.2 Magnitude and surface brightness selection effects 

The DEEP2 spectroscopic survey is limited by apparent R- 
band magnitude of 24.1 AB. This results in a bias against 
redder galaxies at low stellar masses (Willmer et al. 2006). 
However, the stellar masses we use from Bundy et al. (2006) 
are limited by apparent Jf-band magnitude of about 22.5 
AB. This additional K-h&nA selection is close to a selec- 
tion on stellar mass and thus removes much of the bias 
by removing lower-mass blue galaxies from the sample. In 
order to remove any remaining color biases we apply an 
additional redshift dependent lower stellar mass limit of 
log 10 M m i n = 8.5 + z. Our final DEEP2 sample consists of 
~ 800 galaxies with spectro-z's and ~ 2000 galaxies with 
photo-z's. 

A potentially important selection bias for galaxy size 
evolution studies is incompleteness at low surface brightness 
levels, especially at higher redshifts due to the cosmological 



surface brightness dimming, which goes as (1 + z) 4 . To de- 
termine the impact of surface brightness selection effects it 
is customary to create model galaxies, insert them into real 
data, and then run the same source detection algorithms 
that were applied to the real data. For the single orbit 
HST/ ACS F850LP observations from GEMS used by Bar- 
den et al. (2005) and Trujillo et al. (2006), these tests pre- 
dict that surface brightness limits are significant for galaxies 
brighter than the magnitude limit of the COMBO-17 pho- 
tometric redshift survey. However, as shown by Melbourne 
et al. (2007), typical large spiral galaxies at z ~ 1 have 
significant amounts of high-surface-brightness substructure. 
This makes them much easier to detect and measure spec- 
troscopic redshifts for than smooth disks. In this paper 
we use galaxies detected in an effective exposure of two- 
orbits of HST/ ACS imaging (single-orbit each of F606W and 
F814W). Based on the results of Melbourne et al. (2007) we 
conclude that, at the depth of our images, we are not miss- 
ing significant numbers of large star forming galaxies with 
stellar masses greater than ~ 1O 1O M0 at z ~ 1. 

At the other extreme, small galaxies can be confused 
with stars, especially with ground based seeing limited imag- 
ing. However, stars can be reliably separated from galaxies 
based on the fact that stars tend to be brighter than compact 
high-redshift galaxies, and also occupy a different locus from 
galaxies in (B — R, R — I) space (Newman et al. in prep). 
Furthermore, there is nothing to prevent small galaxies from 
being in our photo-z sample. The good agreement between 
the size-mass relations from our spectro-z and photo-z sam- 
ples vindicates the DEEP2 star-galaxy separation method. 
We thus conclude that our sample of star-forming galaxies 
from DEEP2 is unlikely to be affected by biases against very 
small or very large galaxies above the (redshift dependent) 
stellar mass limits of our sample. 

2.2.3 SDSS comparison sample 

As a low redshift comparison sample, we use data from the 
SDSS DR7 (Abazajian et al. 2009). We use a similar method- 
ology to the DEEP2 sample in order to minimize systematic 
biases in the comparison between high and low redshift data. 
Thus, we measure disk half-light radii with 2D bulge+disk 
fits using GIm2d, performed simultaneously on g and r im- 
ages with the bulge and disk sizes constrained to be the 
same in each filter (Simard et al. , in prep). We use galax- 
ies with spectroscopic redshifts from 0.02 < z < 0.20 and 
stellar masses calculated using the relation between r-band 
mass-to-light ratio and (g — r) color from Bell et al. (2003) 
with an offset of -0.1 dex to correspond to a Chabrier (2003) 
IMF. We select blue-cloud galaxies for this study using the 
valley in the (g — r)-stellar mass plane and requiring: 

(g-r)< 0.71 - 0.067(log 10 M stat - 10). (2) 

As with the DEEP2 sample we select galaxies with disk 
light fractions greater than 50% and with disk minor-to- 
major axis ratios greater than 0.5 (corresponding to disk 
inclinations less than 60 degrees for a zero thickness disk). 
We use GIm2d model g and r magnitudes, with k-corrections 
to z = based on SDSS petrosian ugriz magnitudes. 

The SDSS spectroscopic survey is limited by apparent r- 
band magnitude. In the color-stellar mass plane this results 
in a bias against redder galaxies at low stellar masses. In 
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Figure 1. Disk half-light radius-stellar mass (RM) relation of blue-cloud disk-dominated galaxies (bulge fraction less than 0.5) from 
the SDSS. The color corresponds to the observed number density of galaxies on a logarithmic scale. The error bars show the median and 
error on the median in stellar mass bins of width 0.2 dex. The dotted lines show the 84th and 16th percentiles of the size distribution. 
The solid line shows our fit to the median relation using Eq. \3\ For comparison we also show the total half-light radius-stellar mass 
relation of disk-dominated (Sersic n < 2.5) galaxies from Shen et al. (2003, long-dashed line) and the disk half-light radius-stellar mass 
relation from Dutton et al. (2007, short-dashed line). The scatter in the size-mass relation is given in the lower panel. The points show 
our measurements from SDSS, the solid line shows our fit to these data using Eq. [4] while the long-dashed line shows the result from 
Shen et al. (2003). 



order to remove this effect, we apply a minimum stellar mass 
of: log 10 M min = 10.4 + log 10 (z/0.1). Our final SDSS sample 
consists of ~ 130000 galaxies. 

Fig- Q] shows the disk half-light radius-total stellar mass 
relation from our SDSS sample. We fit the median SDSS 
disk size-mass relation with the following double power-law: 

Here a is the slope at low masses (M <C Mo), is the 
slope at high masses (M 3> Mo), and 7 controls the sharp- 
ness of the transition between the two slopes, Mq is the 
transition mass, and Ro is the value of R at Mo. We find 
Iog 10 (Af /[Afg]) = 10.44, log 10 (i? /[kpc]) = 0.72, a = 0.18, 
/3 — 0.52 and 7 = 1.8 provide a good fit to the data. 

We assume the scatter in the SDSS disk size-mass re- 
lation is log-normal with s = <J\ n R\M, with the following 
relation: 

s = s 2 + (si - « a )/[l + (M/Mo) 7 ]. (4) 



Here si is the scatter at low masses, S2 is the scatter at 
high masses, Mo is the transition mass, and 7 controls the 
sharpness of the transition. We find si = 0.47, S2 = 0.27, 
log 10 (M /[M Q ]) = 10.3, and 7 = 2.2 provide a good fit to 
the data. 

For comparison, the total half-light radius-stellar mass 
relation for disk galaxies (Sersic n < 2.5) from the SDSS 
study of Shen et al. (2003) is given by the long-dashed line, 
and the disk half-light size ( 1. 678 x disk scale length)-stellar 
mass relation from Dutton et al. (2007), using the spiral 
galaxy sample of Courteau et al. (2007), is given by the 
short-dashed line. The relation from Dutton et al. (2007) is 
in reasonable agreement with our result considering a sin- 
gle power-law was used by Dutton et al. (2007). The rela- 
tion from Shen et al. (2003) has a shallower slope at high 
masses, and is offset to smaller sizes at all masses. These 
differences can be understood as a consequence of two fac- 
tors: circularization and bulges. The half-light radii used by 
Shen et al. (2003) are measured using circular apertures. 
Since the median major to minor axis ratio of galaxy disks 
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Figure 2. Evolution of the disk half-light radius-stellar mass (RM) relation of disk-dominated blue-cloud galaxies from DEEP2 using 
spectroscopic redshifts. The upper left panel shows the RM relation we derive for disk-dominated blue-cloud galaxies from SDSS. The 
long-dashed line shows the median relation, the dotted lines show Ict scatter. In the other panels the green points show observations 
from DEEP2 using HST imaging in five redshift bins. The median redshift in each bin is given in the upper left corner of each panel, 
together with the central rest-frame wavelength of the imaging used to measure galaxy sizes, and the number of galaxies. The solid red 
line shows a fit to the DEEP2 data by offsetting the z = 0.1 relation (shown by the dashed line) in size. The red circle shows the median 
size for a stellar mass of 3 x 10 10 M o . The offset is given in the upper right corner in each panel, and is summarized in Table [2] The 
red short-dashed lines show a power-law fit to the data in each panel. The vertical dotted lines show the lower stellar mass limit at the 
upper redshift of each redshift bin. 



is ~ 2, a circular size measurement will underestimate the 
true (face-on) half-light radius by a factor of ~ 1.4 (0.15 
dex). This effect explains most of the discrepancy at low 
stellar masses. At high stellar masses there is an additional 
difference, which can plausibly be explained by the increased 
bulge fractions in higher mass galaxies (e.g. Dutton 2009). 
Since bulges tend to be smaller than disks, total half-light 
radii will be smaller than disk half-light radii. 

We find that the scatter in disk sizes is mass dependent, 
with larger scatter at lower masses, in qualitative agreement 
with Shen et al. (2003). At low masses we find that the scat- 
ter of 0.47 in lni? is in agreement with Shen et al. (2003), 
but at high masses we find a scatter of 0.27 in \nR, which 
is smaller than the scatter of 0.34 in In J? reported by Shen 
et al. (2003). As with the difference in the slopes at high 
masses, this difference is plausibly due to use of total half- 
light radii by Shen et al. (2003). We note that in the sim- 
plest ACDM based disk formation models (e.g., Mo, Mao, 
& White 1998) the scatter in disk sizes at fixed mass is 
equal to the scatter in the halo spin parameter, A, which is 
crin A ~ 0.5 (e.g., Bullock et al. 2001b, Macci'o et al. 2007). 
This is in rough agreement with the observed scatter at low 
stellar masses (M s t ar ,$ 10 10 Mq), but at high stellar masses 



(M s tar ^ 10 Mq) the observed scatter is a factor of ~ 2 
lower than predicted by the simple model. This smaller than 
expected scatter has been noted by previous authors (e.g., 
de Jong & Lacey 2000; Shen et al. 2003; Pizagno et al. 2005; 
Dutton et al. 2007). This discrepancy may indicate that mas- 
sive disk-dominated galaxies form in a subset of haloes with 
a biased distribution of halo spin parameters, or that the 
distribution of disk sizes has been modified by secular evo- 
lution (e.g., Shen et al. 2003). 



2.2.4 Evolution 

The evolution of the disk half-light radius-stellar mass rela- 
tion in six redshift bins from z = 0.1 (SDSS) up to redshift 
z ~ 1.2 (DEEP2) is shown in Pigs. H&H The upper left 
panel shows fits to the SDSS disk size-stellar mass relation 
from Fig. □ For both SDSS and DEEP2 the stellar mass is 
the total (i.e., disk plus bulge) stellar mass. The long-dashed 
line shows the median relation, and the dotted lines show 
the la scatter. 

Galaxies at redshifts z > 0.2 from DEEP2 are shown 
as green circles. The solid red lines show fits to the DEEP2 
data obtained by calculating the median offset in size with 
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Figure 3. Evolution of the disk half-light radius-stellar mass (RM) relation of disk-dominated blue-cloud galaxies from DEEP2 using 
photometric rcdshifts. The upper left panel shows the RM relation we derive for disk-dominated blue-cloud galaxies from SDSS. The 
long-dashed line shows the median relation, the dotted lines show la scatter. In the other panels the green points show observations 
from DEEP2 using HST imaging in five redshift bins. The median redshift in each bin is given in the upper left corner of each panel, 
together with the central rest-frame wavelength of the imaging used to measure galaxy sizes, and the number of galaxies. The solid red 
line shows a fit to the DEEP2 data by offsetting the z = 0.1 relation (shown by the dashed line) in size. The red circle shows the median 
size for a stellar mass of 3 X 10 10 M o . The offset is given in the upper right corner in each panel, and is summarized in Table [2] The 
red dashed lines show a power-law fit to the data in each panel. The vertical dotted lines show the lower stellar mass limit at the upper 
redshift of each redshift bin. 



Table 2. Evolution of the disk half-light radius-stellar mass rela- 
tion of blue-cloud disk- dominated galaxies relative to SDSS galax- 
ies at 2 = 0.1, using data from DEEP2 from Figs.[2]&[3] 



redshift range median z A log 10 R50 \ M s tar N 



Spectroscopic Redshifts 



0.2-0.5 


0.38 


-0.05 ± 0.03 


95 


0.5-0.7 


0.61 


-0.08 ± 0.02 


1 15 


0.7-0.9 


0.77 


-0.08 ± 0.03 


267 


0.9-1.1 


0.99 


-0.12 ± 0.02 


192 


1.1-1.4 


1.22 


-0.14 ±0.03 


118 


Photometric Rcdshifts 


0.2-0.5 


0.37 


-0.02 ± 0.01 


360 


0.5-0.7 


0.62 


-0.07 ±0.01 


321 


0.7-0.9 


0.77 


-0.09 ± 0.01 


613 


0.9-1.1 


1.00 


-0.11 ± 0.01 


414 


1.1-1.4 


1.21 


-0.18 ± 0.02 


253 



reference, the 2 = 0.1 relation is shown in each panel with 
black dashed lines. The short-dashed red lines show power- 
law fits to the data in each redshift bin. These agree well at 
M star ~3x 1O 1O M (red dots in Figs. [2]&[3j) but suggest the 
size evolution may be weaker at lower masses and stronger 
at higher masses. 

The evolution of the zero point (size evolution at fixed 
Mstar) is summarized in Fig. [4] The data from DEEP2 
are shown with blue pentagons (open for photo-z, filled for 
spectro-z), and the data from Trujillo et at. (2006) are shown 
with magenta squares. Note that the latter are for total half- 
light radii, not disk half-light radii, but nevertheless the two 
data sets are consistent within the error bars. The DEEP2 
data, however, shows stronger evidence for evolution both 
internally and in comparison to SDSS. For DEEP2 the evo- 
lution of the spectro-z sample is well fit by 



Alog 10 i?50 = 0.018 ±0.002 - 0.44 ± 0.04 log 10 (l±z), (5) 



respect to the z — 0.1 relation. The offsets are given in 

the top right corner of each panel, and in Table [2] For where Alog 10 J?so is the size evolution relative to z = 0.1. 
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Figure 4. Evolution of the disk half-light radius-stellar mass re- 
lation of blue/disk galaxies from redshifts z = 0.1 to z ~ 1.2. 
The pentagons show results for blue-cloud disk-dominated galax- 
ies from DEEP2, where the sizes are the disk half-light sizes. 
The solid blue pentagons show the spectro-2 sample (Fig. 
and the open blue pentagons show the photo-z sample (Fig. [3J. 
Both are well fit by a linear relation (dashed line): A log 10 R50 = 
0.018 — 0.441og 1 g(l + 2). For comparison, the magenta squares 
show data for disk-dominated galaxies (Sersic n < 2.5) from Tru- 
jillo et al. (2006), where the sizes are total galaxy half-light sizes. 
Error bars are lc. 



Figure 5. Comparison of Ho major-axis HWHM radii, i?HWHM> 
of SINS galaxies from Cresci et al. (2009) and Ha half-light radii 
for the same 18 galaxies from Forster-Schreiber et al. (2009). 
The HWHM radii are larger on average by 0.12 dex. Bouche 
et al. (2007) claim that for an exponential distribution, R-hwhm 
equals the exponential scale length, Rd, to within 15%, but then 
iijjwHM should be 0.22 dex smaller than R1/2 (see vertical line 
at left), contrary to what is observed. Thus, either the SINS Ha 
disks are strongly non-exponential and/or the HWHM and half- 
light radii fundamentally disagree. 



2.3 Evolution of the size - rotation velocity 
relation 

The evolution of the size-velocity (RV) relation has been 
studied by Bouche et al. (2007) using data from the SINS 
survey (Forster-Schreiber et al. 2009) at 2 ~ 2 and from 
Courteau (1997) at z = 0. Using half-width half-maximum 
(HWHM) sizes interpreted as exponential disk scale lengths, 
Bouche et al. (2007) found that the RV relation exists at 
2 — 2 with the same zero-point as that at z = 0. 

In this paper we re-determine the zero-point evolution 
using data from SINS at z ~ 2 adding data from Cresci 
et al. (2009) to those of Forster-Schreiber et al. (2009). 
Cresci et al. (2009) determined maximum rotation veloci- 
ties for 18 galaxies at 2 ~ 2 from the SINS survey using 
kinematic modeling of 2D Ha velocity fields. Stellar masses 
were obtained from optical to NIR SED fits assuming So- 
lar metallicity Bruzual & Chariot (2003) stellar population 
synthesis models with a Chabrier (2003) IMF. Galaxy sizes 
were measured from Ha emission maps using two meth- 
ods: Cresci et al. (2009) give HWHM sizes, -Rhwhm, while 
Forster-Schreiber et al. (2009) give circular half-light sizes, 
Ri/2- Both size measurements have been made on the 18 
galaxies from Cresci et al. (2009). 

The HWHM sizes are obtained from a linear Gaussian 
fit to the major axis of the Ha line intensity maps. Ac- 
cording to Bouche et al. (2007), the derived HWHM, once 
corrected for the observed seeing, corresponds to the expo- 



nential scale length, Rd, of the disk. Using simulations of 
model disks these authors claim that Rd measured this way 
is likely to be overestimated by no more than 15%. The 
half-light sizes from Forster-Schreiber et al. (2009) are ob- 
tained from the Ha curves-of-growth measured in circular 
apertures, and corrected for seeing. 

For an exponential disk, Rd = O.6OR1/2, and thus if 
Rhwhm = Rd then Rhwhm should be smaller than -R1/2 
by 0.22 dex, whereas Fig. [S] shows that the opposite is 
true: the HWHM radii are 0.12 dex larger. The fact that 
Rhwhm — 1.3i?i/ 2 was also noted by Forster-Schreiber 
et al. (2009). The discrepancy may be somewhat exagger- 
ated, as the half-light radii are measured through circular- 
apertures, not along the major axis. If the galaxies studied 
by Cresci et al. (2009) were a random sample of galaxies, 
the median disk inclination angle would be 60 degrees, and 
thus the circular radii would be too small on average by a 
factor of ~ 1.4. However, the distribution of disk inclina- 
tions of the galaxies in Cresci et al. (2009) appears skewed 
towards low inclinations: only 3/18 galaxies have inclina- 
tion greater than 60 degrees. The median inclination is 42 
degrees, which implies that the circularization effect is of or- 
der 0.07 dex. In what follows we use both size measurements, 
assuming that -Rhwhm is equivalent to an exponential scale 
length and 1.16 x R1/2 is equivalent to a major-axis half- 
light radius. We will let consistency between various data 
sets determine which, if either, size measurement is more 
likely to be correct. 

Fig. [6] shows the observed evolution of the VMR re- 
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Figure 6. Evolution of the size-stellar mass-rotation velocity relations between redshifts z = and z ~ 2. The ,z = relations are from 
Bell & de Jong (2001) and Dutton et al. (2007), while the z ~ 2 data are from SINS (Cresci et al. 2009; Forster-Schreiber et al. 2009). 
The stellar masses and maximum rotation velocities are from Cresci et al. (2009). In the upper panels the size is the exponential disk 
scale length (from Cresci et al. 2009), whereas in the lower panels it is the half-light radius (from Forster-Schreiber et al. 2009). These 
two size measurements give different, and inconsistent, results for the evolution of disk-galaxy sizes. 



lations from z ~ 2 to z ~ 0. The SINS data at z ~ 2 
are shown with circles. Two versions of radii are shown: 
HWHM Ha radii by Cresci et al. (2009) (upper panels) and 
(de-circularized) half-light Ha radii by Forster-Schreiber 
et al. (2009) (lower panels). The redshift z — relations 
are shown as dashed lines and are obtained as follows. We 
use the M sta r — V m ax relation from Bell & de Jong (2001) 
derived from /\-band luminosities using a mass dependent 
extinction correction, subtracting 0.1 dex from the stellar 
masses to correspond to a Chabrier (2003) IMF, and using 
a bisector fit: 



log'i 



= -0.61 + 4.51 lc^ 



V m 



[100 kms- 1 ] 



(6) 



We note that the velocity used by Bell & de Jong (2001) 
is not Knax, but actually Vkat, the velocity in the flat outer 
part of the rotation curve. However, as shown by Dutton 
et al. (2010b) the TF relation derived using V^ a x for these 
galaxies is identical to the one derived using Vfiat- 

We use the local /-band disk scale length-stellar mass 
relation from Dutton et al. (2007) which used an orthogonal 
fit: 



log -ra = a348 + - 2811o ^Iic^ 



(7) 



where we have adopted stellar masses from Bell et al. (2003) 
-0.1 dex. Taking these MV and RM relations, we infer the 
mean relation between Rd,i and Vm^: 



log'i 



Rd 



[kpc] 



0.177 + 1.267 log 1 



V m 



[100 kms" 1 ] 



(8) 



In order to compare with disk half-light radii we use R$o ,i = 
1.678-Rdj- For simplicity we use the /-band sizes to compare 
with the z ~ 2 data. But we note that the evolution would 
not change significantly, in fact it would be slightly stronger, 
if we used the V-band RM relation that we derived from the 
SDSS. 

Due to the small numbers of galaxies in the z ~ 2 sam- 
ple we do not attempt to fit the slopes of the high-redshift 
relations. Instead we fix the slopes to the z — values and 
measure the median offset of the z ~ 2 galaxies with respect 
to the z — relations. The values of these offsets are given 
at the bottom of each panel. 

The upper right panel shows the disk scale length- 
maximum rotation velocity relation using the Cresci 
et al. (2009) SINS data. The z ~ 2 data have a marginally 
higher zero point (0.09 + 0.05 dex) compared to the redshift 
zero data, which is in agreement with the findings of Bouche 
et al. (2007). Given that the TF relation evolves to lower 
stellar masses (—0.44 + 0.08 dexj^at fixed velocity at higher 
redshifts (lower left panel), this implies that at fixed stellar 



4 Note that the difference in evolution of 0.05 dex between our 
result and that of Cresci et al. (2009), who use the same data, is 
caused by our different normalization of the 2 = relation from 
Bell & de Jong (2001). In order to convert the stellar masses 
derived assuming a "diet"-Salpeter IMF of Bell & de Jong (2001) 
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mass, sizes should be larger at higher redshifts. As is shown 
in the upper middle panel, this is indeed the case, with sizes 
larger by 0.20 ± 0.05 dex at z ~ 2 than z — 0. The evolution 
of the R d - M star relation from the SINS Cresci et al. (2009) 
data is thus of opposite sign to that of previous studies (see 
§ 12. 2[) . In particular it is inconsistent with the z ~ 2 results 
from Trujillo et al. (2006) and Williams et al. (2010), who 
found that disk/star forming galaxies galaxies were a factor 
of ~ 2 smaller at fixed stellar mass. 

The lower middle and right panels show the RM and 
RV relations using the SINS curve-of-growth Ha half-light 
radii from Forster-Schreiber et al. (2009). These panels show 
that the half-light sizes are smaller at a given stellar or max- 
imum rotation velocity, relative to z = galaxies. The offset 
for the RM relation is -0.07 ± 0.05 dex, while for the RV 
relation it is -0.22 ± 0.06 dex. The evolution in the RM 
relation is still however, weaker than obtained by Trujillo 
et al. (2006), who find —0.28 dex. As we discuss below in 
§ 13.31 this remaining difference can be accounted for by the 
difference between half-light radii measured in Ha and those 
measured in rest-frame /-band light. 



3 THEORETICAL EXPECTATIONS 

We now compare these observational results to the predic- 
tions of various theoretical models. 



3.1 Dark matter haloes 

The zeroth order prediction for the evolution of the VMR 
scaling relations of disk galaxies is given by the evolution of 
the VMR relations of dark matter haloes. This is obtained 
by assuming that: 

(i) The total mass profile of the galaxy and halo is isother- 
mal; 

(ii) The galaxy mass fraction, m ga i = M ga i/M v i r , is a con- 
stant and independent of redshift. Here M ga i is the galaxy 
mass (i.e., the sum of stellar mass and cold gas), and M v i r 
is the total mass within the halo virial radius; 

(in) The galaxy spin parameter, A ga i = (j ga i/m ga i)A, is 
a constant and independent of redshift (for a given halo). 
Here j ga i = J ga i / Jvir is angular momentum fraction, where 
J ga i is the total angular momentum of the galaxy and J v i r 
is the total angular momentum within the virial radius, and 
the halo spin parameter, A, is given by 



A 



Jvir I SI 1/2 



GM 



5/2 



Jvir/Mvir y 1 /^ 
V2 -RvirVvir 



(9) 



Here 7? v ir is the virial radius, V V ir is the circular velocity at 
the virial radius, E is the halo's energy, and f c = 1. For 
more general haloes f c measures the deviation of E from 
that of a singular isothermal sphere. For an NFW halo f c ~ 
| + (c/21.5) ' 7 (see Mo, Mao, & White 1998), where c is the 
halo concentration; 

(iv) The galaxy is 100% stars, i.e., there is no cold gas. 



into a Chabrier (2003) IMF, which is used at z ~ 2, we adopt -0.1 
dex, whereas Cresci et al. (2009) adopted -0.15 dex. 



We refer to this model as the singular isothermal sphere 
(SIS) model. Under these assumptions the disk size scales as 
the size and circular velocity of the halo (e.g., Mo, Mao, & 
White 1998): fla.sis — XRvir/V^ x V v i r , and the stellar mass 
scales with the halo mass: M star oc M v j r . 

Assuming the halo is defined by spherical top hat col- 
lapse: Mvir = (4/3)7ri?J ir AvirPcrit, then the evolution of the 
MV, RM, and RV relations is given by: 



Mvir (z) 


Vvir(z)" 


3 


'Avir(z)" 


-1/2 


\H(z)] 


Mvir.O 


Vvir.O 




A v ir,0 




H 



Rvir{z) _ 


'M vil (z)~ 


1/3 


~A vir (2)~ 


-1/3 


\H(z)] 


Rvir,0 


M v i r ^o 




A v ir,0 




Ho 



gvir(g) 
-Rvir.O 



-2/3 



Vvh(z) 



Vvir.O 



Ayjr(.z) 
Avir.O 



-1/2 



H(z 



Ho 



(10) 

(11) 
(12) 



The zero point evolutions are obtained by setting the leading 
terms on the RHS=1. Thus the evolution is governed by 
the evolution of the Hubble parameter H(z) and the halo 
overdensity within the virial radius, A v i r (z). The evolution 
of the Hubble parameter (in a flat cosmology) is given by 



H(z) =# [Oa+«m(1 + z) ] 



3U/2 



(13) 



For the evolution of the halo overdensity, A v - a (z), we use 
the fitting formula from Bryan & Norman (1998): A v i r = 
187T 2 + 82x — 39a; 2 , where x = Q,(z) — 1 and Q(z) is defined 
asfi M (l + 2) 3 [H(z)/H ]- 2 . 

Thus to first order the evolutions are determined by 
the evolution of the Hubble parameter and to second or- 
der by the evolution of the halo overdensity. These pre- 
dicted evolutions, for a ACDM cosmology with Qm = 0.3 
and Oa = 0.7, are shown by the dashed lines in Fig. [7] 
For the M v i r — Vvir and i? v ir — Vvir relations the evolution 
scales roughly as (1 + z) -1,3 , while for the R v i T — M V - 1T re- 
lation the evolution scales roughly as (1 + ,z)~ ' 8 . The ob- 
served evolution of the maximum circular velocity - stel- 
lar mass - half-light radius relations is given by the col- 
ored symbols. For the stellar mass-velocity relation we use 
data from: Conselice et al. (2005, green triangles); DEEP2 
- Kassin et al. (2007, blue pentagons); and SINS - Cresci 
et al. (2009, red open circle). For the size-stellar mass re- 
lation we use data from: Trujillo et al. (2006, magenta 
squares); SINS - Cresci et al. (2009, red open circle) and 
Forster-Schreiber et al. (2009, red open square); Williams 
et al. (2010, green triangles); and DEEP2 - this paper (blue 
pentagons). For the size- velocity relation we use data from: 
SINS - Cresci et al. (2009, red open circle) and Forster- 
Schreiber et al. (2009, red open square). These observations 
are described in more detail in Sj2] As has been noted by 
previous authors, the observed evolution of all three VMR 
relations is weaker than this simple prediction from dark 
matter haloes. 

The SIS model makes a number of simplifying assump- 
tions that are likely to be incorrect. Firstly, total mass den- 
sity profiles are not expected to be globally isothermal. Dark 
matter haloes in cosmological simulations have V max > V v - U 
(e.g., Bullock et al. 2001a). In addition, the contribution of 
the baryons to the inner potential can increase the observed 
V m ax even further. 

Using the models of Mo, Mao, & White (1998), 
Somerville et al. (2008) showed that including the expected 
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Figure 7. Evolution of the circular velocity - mass - size relations in observations and models. The data are given by colored symbols, 
as indicated, and are described in J2] The evolution of the stellar relations (V m ax — Afstar — Rsov) m a ACDM disk galaxy evolution 
model are given by solid lines. The models and data are consistent, with the exception of the sizes of galaxies in the SINS survey as 
measured by Cresci et al. (2009), which appear to be a factor of ~ 3 too large (red open circles). This discrepancy is partially removed 
when using Ha half-light sizes from Forster-Schreiber et al. (2009) (red open squares), and fully removed after correcting these sizes 
from Ha to rest- frame /-band (black solid circles). For comparison, the evolution of the baryonic (V^nax — Af ga i — i? ga i) and halo virial 
(Vvir — Mvii — Rvir) relations are given by the dotted and dashed lines, respectively. The evolution of the theoretical virial, baryonic, 
and stellar relations are given in Table [3] 



evolution in dark halo concentrations, as well as the contri- 
bution of baryons to the potential, results in weaker evolu- 
tion than predicted by the SIS model. However, Somerville 
et al. (2008) make the assumption that disks have no gas, 
and thus, that the stellar mass is equal to the baryonic mass 
and the V-band half-light radius is equal to the baryonic 
half-mass radius. In order to determine whether the stellar 
relations evolve differently from the baryonic relations, one 
needs to follow the evolution of stellar and gas disks in a 
self-consistent way. This requires realistic cosmological sim- 
ulations with gas or semi-analytic models (SAMs). Here we 
use the latter. 



3.2 Disk galaxy semi-analytic model 

In order to calculate the evolution of stellar disks as opposed 
to baryonic disks, we use the disk galaxy SAM of Dutton 
& van den Bosch (2009). This model consists of disks that 
grow inside evolving NFW haloes (Navarro, Frenk, & White 
1997) , with structure determined from cosmological N-body 
simulations (Bullock et al. 2001a; Maccio, Dutton & van 
den Bosch 2008), a median spin parameter A, which is in- 
dependent of redshift, and halo specific angular momentum 
distributions from Sharma & Steinmetz (2005), which are 
specified by a parameter, a. The cumulative distribution of 
specific angular momentum P(< s), where s is the specific 
angular momentum in units of the total specific angular mo- 
mentum, is given by 

P(<s) = y(a,as), (14) 

where 7 is the incomplete gamma function. In this model, 
disks are not formally exponential, but the stellar disks can 
often be well described by an exponential profile over several 
scale lengths (Dutton 2009). 

To build a model as close as possible to the SIS model, 
we maintain assumptions (ii) and (iii) from above. As in 



Somerville et al. (2008), the baryonic disk is in dynamical 
equilibrium inside an NFW halo which evolves with redshift 
according to cosmological simulations (Bullock et al. 2001a). 
The key difference is that we do not make assumption (iv), 
i.e., that the disk is 100% stars. In our model, the evolu- 
tion of the stellar and gas disks (and hence stellar and gas 
mass) is governed by the radial variation of star formation, 
gas recycling, and accretion. The stellar mass is thus always 
less than the baryonic mass, and the stellar disk is usually 
smaller in radius than the baryonic disk. Because our model 
follows stellar populations as a function of radius, we calcu- 
late the sizes of model galaxies in optical to NIR light, as well 
as stellar and baryonic mass. This allows us to test whether 
the sizes measured in different pass-bands and masses are 
equivalent. 

An additional difference between our models and those 
of Somerville et al. (2008) is halo contraction. Somerville 
et al. (2008) assumed that haloes contract according to the 
Blumenthal et al. (1986) adiabatic contraction model. In our 
model, we leave the haloes uncontracted, as models with 
dark halo contraction (and standard IMFs) are unable to 
reproduce the zero points of the VMR relations (Dutton 
et al. 2007). There are a variety of astrophysical processes 
that could reverse the expected effect of halo contraction. 
These include dynamical friction from massive clumps of 
baryons (e.g., El-Zant, Shlosman & Hoffman 2001; Mo & 
Mao 2004; Elmegreen et al. 2008; Jardel & Sellwood 2009), 
dynamical friction due to bars (e.g., Weinberg & Katz 2002; 
Sellwood 2008), and rapid mass outflows due to supernova 
feedback (e.g., Mo & Mao 2004; Read & Gilmore 2005; Gov- 
ernato et al. 2010), so our choice of not including adiabatic 
contraction is at least physically plausible as well as being 
empirically motivated. As we will see from the success of 
our models in matching the evolution of the VMR relations, 
halo contraction is also not required in order to reproduce 
the observed evolution of the VMR relations. 
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Table 3. Evolution of the stellar, baryonic, and virial velocity-mass-size relations of our 
disk-galaxy evolution model, with m ga i = 0.04 and A = 0.035, relative to redshift z = 
as shown in Fig. [7] 



2 


login (1 + *) 


A log 10 Mstar 


Vmax 


Alog 10 Mg al |Vmax 


Alog 10 Af vir |y vir 


0.00 


0.000 


0.000 




0.000 


0.000 


0.10 


0.041 


-0.001 




0.006 


-0.035 


0.30 


0.114 


-0.027 




-0.014 


-0.106 


0.51 


0.179 


-0.065 




-0.039 


-0.177 


0.70 


0.230 


-0.105 




-0.066 


-0.238 


1.01 


0.303 


-0.171 




-0.115 


-0.327 


1.41 


0.382 


-0.264 




-0.190 


-0.435 


2.00 


0.477 


-0.418 




-0.307 


-0.568 


2.50 


0.544 


-0.543 




-0.394 


-0.665 


3.00 


0.602 


-0.661 




-0.472 


-0.749 


4.00 


0.699 


-0.872 




-0.610 


-0.892 


2 


iogioC 1 + z ) 


A log 10 R 50V 


Af s tar 


A logi #50gal|AW 


Alog 10 i? vir |M vir 


0.00 


0.000 


0.000 




0.000 


0.000 


0.10 


0.041 


0.003 




0.012 


-0.024 


0.30 


0.114 


-0.019 




-0.001 


-0.071 


0.51 


0.179 


-0.046 




-0.022 


-0.118 


0.70 


0.230 


-0.071 




-0.042 


-0.159 


1.01 


0.303 


-0.116 




-0.071 


-0.218 


1.41 


0.382 


-0.175 




-0.112 


-0.290 


2.00 


0.477 


-0.270 




-0.170 


-0.379 


2.50 


0.544 


-0.339 




-0.224 


-0.443 


3.00 


0.602 


-0.402 




-0.276 


-0.500 


4.00 


0.699 


-0.519 




-0.368 


-0.595 


2 


logio(l + *) 


A logi r mv 


Vmax 


A log 10 i?5n ga l|Vmax 


Alog 10 i? vir |V v ir 


0.00 


0.000 


0.000 




0.000 


0.000 


0.10 


0.041 


0.003 




0.010 


-0.035 


0.30 


0.114 


-0.035 




-0.016 


-0.106 


0.51 


0.179 


-0.077 




-0.046 


-0.177 


0.70 


0.230 


-0.120 




-0.069 


-0.238 


1.01 


0.303 


-0.192 




-0.121 


-0.327 


1.41 


0.382 


-0.291 




-0.197 


-0.435 


2.00 


0.477 


-0.448 




-0.306 


-0.568 


2.50 


0.544 


-0.557 




-0.395 


-0.665 


3.00 


0.602 


-0.672 




-0.479 


-0.749 


4.00 


0.699 


-0.850 




-0.614 


-0.892 



In the model used here, we adopt a cosmology with 
(n M ,{lA,h,a$,n) = (0.3,0.7,0.7,0.8,1.0), which is close to 
that of the WMAP 5th year results (Dunkley et al. 2009). 
We adopt a galaxy mass fraction of m ga i = (Af s t ar + 
Mgas)/M v i r = 0.04, a median spin parameter of A ga i = A = 
0.035, and a median angular momentum shape parameter 
of a — 0.9. These parameters are motivated by observa- 
tions of m ga i (Hoekstra et al. 2005; Dutton et al. 2010b) and 
theoretical predictions for A (Bullock et al. 2001b; Maccio 
et al. 2007) and a (Sharma & Steinmetz 2005). They also 
result in models that roughly reproduce the observed VMR 
relations at z = (Fig. [8J. However, this model does not 
exactly reproduce the slopes of the local VMR relations. 
Doing so requires either m ga i or A ga i to vary with halo mass 
(e.g., Shen et al. 2003; Dutton et al. 2007). 

We generate a Monte Carlo sample of 2000 galaxies 
with halo masses between M vil = 10 10 ' 3 - W 13 5 h~ 1 M Q , 
log-normal scatter in spin parameter of a\ n \ =0.5 (Bul- 
lock et al. 2001b), log-normal scatter in halo angular mo- 



mentum profile of <Ti na = 0.25 (Sharma & Steinmetz 2005), 
and log-normal scatter in halo concentration of ai nc = 0.25 
(Maccio et al. 2008). We determine the evolution of the zero 
points of the model VMR relations by fitting the MV and 
RV relations for 2.1 ^ log 10 (V max /[kms -1 ]) ^ 2.5 and the 
RM relations for 9.0 < log 10 (M star /[A/ Q ]) < 11.0. We then 
calculate the evolution at log 10 (V max /[ kms -1 ]) = 2.3 and 
log 10 (M at „/[M©]) = 10.5. 

The solid lines in Fig.[7]show evolution in the zero points 
of the Vmax — Mstar — R&ov (i.e., maximum circular veloc- 
ity, stellar mass, V-band half-light radius) relations of this 
model from redshifts z — 4 to z = 0. These relations show 
weaker evolution than the virial relations (dashed lines), but 
stronger evolution than the baryonic relations (dotted lines) . 
The evolution of the theoretical virial, baryonic, and stellar 
relations are given in Table [3] 

One interpretation of the weak evolution of the galaxy 
scaling relations since z ~ 1 is that individual galaxies also 
evolve weakly since z ~ 1. However, disk-galaxies were form- 
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Figure 8. Evolution of the maximum circular velocity - stellar mass - half-light radius relations in our ACDM disk-galaxy evolution 
model. Power-law fits to the (V max — Mstar ~ -Rgol) relations at redshifts 2 = 4, 2, 1, are given by solid black lines. The 2 = observations 
are given by short-dashed black lines. The symbols show the evolution of an individual galaxy with redshift 2 = halo virial mass of 
M v ; r = lO 12 h _1 M0 at 2 = 4, 3, 2, 1, (from lower left to top right). The black circles show the evolution of V max — Af s t ar — Rsoi- There 
is substantial growth in these three quantities, in a roughly power- law from 2 = 4 to 2 = (black dotted lines). The green triangles 
show the evolution of the SIS model, which is identical to the evolution of the halo virial quantities. To facilitate comparison between 
the galaxy and virial properties, the red squares show M s t ar — V vlr and M s t ar — R^ gig. The weak evolution of galaxy scaling relations 
since 2 ~ 1 is due to the fact that individual galaxies grow roughly along the scaling relations. 



ing stars at much higher rates in the past (e.g., Noeske 
et al. 2007a), and thus stellar masses are expected to grow 
significantly since z ~ 1. Simple models for the evolution 
of the star formation rate-stellar mass relation suggest that 
galaxies with present day stellar masses of 3 x 10 10 Mq had 
stellar masses a factor of ~ 2.5 lower at 2 = 1 (Noeske 
et al. 2007b). This amount of evolution in stellar masses is 
consistent with that predicted by our model. Fig. [8] shows 
the evolution of a galaxy with present day stellar mass of 
3 x 10 10 M Q . Since z = 1 its stellar mass has increased by 
a factor of ~ 2.5, its half-light size has increased by a fac- 
tor of ~ 2, and its maximum circular velocity has increased 
by a factor of ~ 1.15. In terms of samples of galaxies, the 
evolution since z = 1 is just a factor of ~ 1.5 increase 
in stellar masses, at fixed Vm ax , and a factor of ~ 1.3 in- 
crease disk size, at fixed stellar mass. Thus, the scaling re- 
lations between properties of galaxies (V max — M star — R501) 
evolve only weakly since z ~ 1 because individual galaxies 
evolve roughly along the scaling relations. Similar theoret- 
ically based conclusions have been made previously for the 
stellar mass - velocity relation (Portinari & Sommer-Larsen 
2007) and the size-stellar mass relation (Firmani & Avila- 
Reese 2009). 



The differences between the evolution of the baryonic 
and virial MV relations (dotted vs. dashed lines in Fig. 0) 
are due solely to evolution in the ratio between Vm ax and Vvir 
and not due to an evolution in the ratio between the bary- 
onic mass and virial mass, m ga i, which is fixed to a constant 
in this model. The ratio between Vm ax and V v i T increases 
towards lower redshift due to higher concentrations in lower 
redshift haloes. This results in less evolution in V max than 
V v ir at fixed M ga i or M v i r , and hence less evolution in M ga i 
or Mvir at fixed Vm ax than at fixed Vvir- The evolution of 
Kn ax and Vvir for an individual galaxy with redshift 2 = 
virial mass of M v i r = 10 12 /i _1 Mq and median halo param- 
eters is shown in the left panels of Fig. [8] This shows that 
for high redshifts (z > 2) Vmax — Vvir, but for low redshifts 
(z < 1) Vvir remains constant while V max continues to in- 
crease. At redshift 2 = 0, V max — 1.3Vvir, which is consistent 
with recent measurements (Dutton et al. 2010b). The lower 
middle panel of Fig. [8] shows that i? v ir and M v ir both con- 
tinue to increase at low redshifts. The reason V v i r remains 
roughly constant while _R v ir and M v i r increase is due to a 
trade off between the addition of new mass (which increases 
Vvir) and the increase of the virial radius (which decreases 
Vvir because the circular velocity of the NFW profile declines 
at large radius). 
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Since the baryonic mass fraction we adopt is only ~ 25% 
of the universal baryon fraction, it is certainly plausible 
that the mechanisms responsible for making galaxy forma- 
tion inefficient result in m ga i varying with redshift. How- 
ever, to first order, variation in m ga i moves galaxies along 
the VM relation (e.g., Navarro & Steinmetz 2000; Dut- 
ton et al. 2007), and thus we expect that large changes in 
m ga i with redshift would be needed in order to significantly 
change the evolution from that predicted by our constant 
m ga i model. 

The differences between the evolution of the baryonic 
and virial RM relations in Fig. [7] are caused by evolu- 
tion in the ratio between baryonic and virial sizes (be- 
cause the baryon mass fraction is a constant), which are 
roughly given by: -Rso ga i /Rvir oc A ga l(Kir/K ia x)/c 1/2 ■ Thus 
as with the MV relations, the differences between the bary- 
onic and virial RM relations are driven by the evolution in 
V^iax/Vvir. Unlike the VM relation, the baryonic RM rela- 
tion is sensitive to the adopted baryon mass fraction (Dut- 
ton et al. 2007). Thus if m ga i decreases (or increases) with 
increasing redshift, this will result in weaker (or stronger) 
evolution in the baryonic RM relation. 

The differences between the evolution of the baryonic 
and virial RV relations are determined by evolution in both 
Rsogai/Rvir and Vm&x/Vvir- Evolution in the size ratio ac- 
counts of ~ 2/3 of the evolution difference, while evolution 
in the velocity ratio accounts for the remaining ~ 1/3. As 
with the RM relation, the RV relation is sensitive to the 
galaxy mass fraction, with a similar sign dependence. 

The differences between the stellar and baryonic MV 
relations in Fig.[7|can be understood as a result of the higher 
cold gas fractions at higher redshifts in our model. We note 
that the evolution in the cold gas fractions in our model is 
relatively modest, with a factor of < 2 increase in gas masses 
at fixed stellar mass between redshift z = and z — 2. As 
shown by Dutton, van den Bosch & Dekel (2010a), more 
general models (including cooling and outflows) also predict 
weak evolution in cold gas fractions, but strong evolution in 
molecular gas fractions. Both of these predictions are con- 
sistent with recent observations (e.g., Erb et al. 2006; Daddi 
et al. 2010; Tacconi et al. 2010; Puech et al. 2010). 

Higher cold gas fractions also contribute to different 
evolution in the stellar and baryonic RV and RM relations. 
But there are also contributions from differences between 
sizes in stellar mass and sizes in optical light, which are a 
consequence of the inside-out nature of stellar disk growth 
in our models. These differences are discussed in more detail 
below. 



3.3 Reconciling SINS sizes with other 
observations 

Our model nicely reproduces the observed zero point evo- 
lution of the stellar mass Tully-Fisher (M star — Kn ax ) and 
size-stellar mass relations from z — 2.2 to z = (Fig. [7]). The 
model also predicts a factor of ~ 3 evolution in the zero point 
of the optical half-light size-velocity relation (Rsov — Vmax) 
over this redshift range (solid black line in the right panel of 
Fig. [7]). However, the Rd — Vm ax data from SINS (using disk 
scale lengths from Cresci et al. 2009) indicate weak evolution 
in the opposite direction from z = 2.2 to z — (red open cir- 



cle in the right panel of Fig. [7]). Recall that we have derived 
the evolution of the Rd — Vin ax relation by comparing with 
the /-band Rd — Knax relation at z = (see § [2T3J) . Since the 
Cresci et al. (2009) data set also implies that higher redshift 
galaxies are larger at fixed stellar mass (red open circle in the 
middle panel of Fig. [7|, this obviously means that the evo- 
lution of the scaling relations between different data sets are 
inconsistent. The offset between the observed RM evolution 
from Trujillo et al. (2006) and the Cresci et al. (2009) SINS 
data is also a factor of ~ 3. Thus if the Cresci et al. (2009) 
disk scale lengths of the SINS galaxies could be reduced by 
a factor of ~ 3, then all of the data sets would be consis- 
tent. Furthermore, the data would be consistent with our 
simple, but cosmologically motivated, model for disk-galaxy 
evolution. 

As discussed in §2.31 there are two size measurements 
now available for SINS galaxies, both based on Ha imag- 
ing: major axis HWHM sizes by Cresci et al. (2009) (which 
have been interpreted as being equivalent to disk scale 
lengths), and circular half-light radii by Forster-Schreiber 
et al. (2009). In flOlwe calculated the evolution of the SINS 
data at z ~ 2.2 relative to the J-band RM and RV relations 
at z = 0. We found that the two size measurements resulted 
in very different amounts of evolution, with the half-light 
sizes giving a factor of ~ 2 smaller sizes at higher redshift. 
Thus using the SINS half-light sizes by Forster-Schreiber 
et al. (2009) (red open squares in Fig.0 removes most of the 
discrepancy that exists when using the Cresci et al. (2009) 
HWHM sizes interpreted as disk scale lengths. 

However, there is still a small discrepancy at the 0.2 dex 
level between the evolution of the SINS half-light sizes with 
those from Trujillo et al. (2006), Williams et al. (2010), and 
our models. Below we discuss whether this difference can be 
explained by the difference between sizes measured in Ha 
compared to rest-frame /-band light. Since the SINS survey 
is not a volume limited sample of galaxies at high-redshift, 
and at least some of the SINS sample was selected on the 
basis of disky morphology with spatially resolved velocity 
gradients, another possibility is that a selection bias against 
small galaxies exists in the SINS survey. 

3.4 The relation between sizes in Ha vs. 
rest-frame optical light 

How well do galaxy sizes measured in Ha trace those mea- 
sured in rest-frame optical stellar light? We now address this 
question using the disk SAM discussed in § 13.21 To make our 
comparison, we assume that Ha is a reliable tracer of recent 
star formation, where recent star formation in our models 
is the star formation within the last time step (i.e., ~ 30 to 
100 Myr). For simplicity, we also ignore the effects of dust, 
which may modify the relation between sizes measured in Ha 
and optical light. Fig.[9]shows the results where we compare 
the half- mass radii of recent star formation, Rsf, in our disk 
formation model, to the half-light radii in the near-UV to 
near-IR pass bands. 

We find that the half-mass radius of the recent star 
formation is on average a factor of ~ 2 larger than the half- 
mass radius of the total stellar mass. This is a signature of 
the "inside-out" nature of stellar disk growth in our models. 
The inside-out growth is due to a combination of two fac- 
tors: (1) For an individual galaxy the baryonic disk grows 
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Figure 9. Ratios of galaxy sizes (half-mass or half-light) in recently formed stars (within the last lOOMyr), Rsf, to sizes in various 
rest-frame optical-near-infrarcd pass bands, Rx, at redshift 2 = 2 (left) and as a function of redshift (right) in a cosmological disk galaxy 
evolution model. The numbers show the median ratio together with the 16th and 84th percentiles of the distribution. The half-mass 
radius of the recently formed stars is, on average, a factor of ~ 2 larger than the half-mass radius of the total stellar mass. This reflects 
the fact that stellar disks are growing inside out in this model. These differences persist, but are not as big, when comparing i?SF to 
optical half-light radii. At z = 2 the rest-frame /-band half-light sizes are ~ 0.15 dex smaller than star formation sizes. These ratios 
evolve only weakly with redshift in our model. 



with time due to an increase in the specific angular momen- 
tum. This is cosmological inside-out growth; (2) The density 
dependence of the star formation law (star formation is less 
efficient at lower gas densities). This is star formation in- 
duced inside-out growth. 

Using cosmological hydrodynamical simulations of 
galaxies at redshift z ~ 2 Sales et al. (2009) found a similar 
factor of ~ 1.8 difference between the half-mass sizes of the 
dense (star forming) gas and the half-mass sizes of the stellar 
mass. However, as we go from stellar mass sizes to NIR sizes 
to optical sizes, the differences with respect to the star for- 
mation sizes decrease. In terms of optical half-light radii, the 
star formation size is just 0.1 dex higher than V-band sizes 
and 0.15 dex higher than /-band sizes at z = 2. Thus Fig. [9] 
suggests that the differences between Ha and optical sizes 
by themselves are unlikely to explain the full factor of 3 dis- 
crepancy between the SINS disk scale lengths and the other 
observations shown in Fig. [7] But, as shown by the solid red 
circles, these are of the right magnitude to fully remove the 
lingering small discrepancy between the Forster-Schreiber 
et al. (2009) half-light radii and the Trujillo et al. (2006) 
data and our models. 

To summarize, we have shown that the radii of SINS 
galaxies can plausibly be reconciled with other observational 
data and with theoretical models provided proper radii are 
used (half-light radii, not HWHM) and provided a small cor- 
rection is applied to convert Ha radii to optical radii. Ideally 
this agreement should be verified directly with sizes mea- 
sured in rest-frame optical light, and preferably using the 
same techniques and definitions as used at lower redshifts. 
It would also be desirable to test this concordance by mea- 



suring the evolution of the VMR relations at redshifts z > 1 
using larger samples of galaxies and over a wider range of 
masses than current studies. The measurement of rest-frame 
optical sizes of galaxies at z ~ 2 is currently possible with 
//-band imaging with HST, and H & K-band imaging with 
adaptive optics from the ground. 

The right panel of Fig. [9] predicts that not much evo- 
lution is expected in the size ratios between different pass- 
bands. Thus as a test of how realistic the size ratios of our 
models are at high redshift, we can compare them to ob- 
servations of size ratios at low redshift. Fig. [9] shows that 
at redshift z = the sizes in the B-band should be larger 
than those in the V, R and //-bands by 0.037, 0.061, and 
0.145 dex, respectively. This is consistent with observations 
of nearby spiral galaxies by MacArthur, Courteau & Holtz- 
man (2003) who find differences of 0.029, 0.061, and 0.127 
dex. While this does not prove that our model predicts the 
correct size ratios for z — 2 galaxies, it does provide indirect 
support. 

Another prediction of the model is that sizes in molec- 
ular gas should be only slightly larger than sizes in Ha or 
UV light. This is a consequence of the almost linear rela- 
tion between star formation rate surface density and molec- 
ular gas surface density in our model (See Fig. Al, Dutton 
et al. 2010a). At higher molecular gas densities the slope of 
the star formation law in our model is steeper than unity, 
with an asymptotic value of 1.4, in agreement with the stan- 
dard Kennicutt-Schmidt relation (Kennicutt 1998). Since 
molecular gas disks are predicted to be higher density at 
higher redshifts, this results in a slight evolution in the ratio 
between molecular gas sizes and star formation rate sizes 
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Table 4. Evolution of the TF relation in optical to near-IR luminosities from disk-galaxy evolution model relative to redshift z = as 
shown in Fig. 1101 
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0.00 


0.000 


0.000 


0.000 


0.000 


0.000 


0.000 


0.10 


0.041 


0.123 


0.068 


0.022 


0.001 


-0.005 


0.30 


0.114 


0.236 


0.151 


0.076 


0.034 


0.021 


0.51 


0.179 


0.296 


0.200 


0.107 


0.051 


0.034 


0.70 


0.230 


0.329 


0.225 


0.123 


0.060 


0.038 


1.01 


0.303 


0.357 


0.250 


0.137 


0.060 


0.036 


1.41 


0.382 


0.374 


0.263 


0.139 


0.048 


0.018 


2.00 


0.477 


0.353 


0.234 


0.094 


-0.010 


-0.051 


2.50 


0.544 


0.319 


0.191 


0.043 


-0.070 


-0.118 


3.00 


0.602 


0.279 


0.143 


-0.016 


-0.132 


-0.186 


4.00 


0.699 


0.192 


0.047 


-0.123 


-0.245 


-0.311 



in our model (right panel, Fig. [9]). The similarity between 
molecular gas sizes and UV sizes has already been observed 
for a handful of star forming galaxies at redshift z ~ 1.5 
(Daddi et al. 2010), which provides further support for our 
model. 

3.5 The evolution of the Tully-Fisher relation in 
optical to near-IR luminosities 

The evolution of our model Tully-Fisher relations in opti- 
cal to near-IR luminosities is shown in Fig. [10] and in Ta- 
ble SI Contrary to the stellar and baryonic TF relations, 
which show a decrease in mass at fixed circular velocity (solid 
and dotted lines), the luminosity TF relations show an in- 
crease in luminosity at fixed circular velocity between today 
and redshift z ~ 1. Luminosities in bluer pass-bands show 
stronger evolution than redder pass-bands. In the B-band, 
the evolution is 0.36 dex (or 0.9 magnitudes) since 2=1, 
while in the A'-band, the evolution is just 0.04 dex (or 0.1 
magnitudes) since z = 1. A similar theoretical prediction 
of 0.85 mag for the evolution of the B-band TF relation 
since 2 = 1 has been shown by Portinari & Sommer-Larsen 
(2007) using cosmological hydrodynamical simulations. Us- 
ing a semi-analytic model Firmani & Avila- Reese (2003) find 
differences between the evolution of the B-band and H-band 
TF relations which are qualitatively similar to what we find 
in our models, though in detail there are differences in the 
absolute evolution. 

Most observations find no or weak evolution in the J- 
and A"-band TF relations since z ~ 1 (Conselice et al. 2005; 
Weiner et al. 2006b; Flores et al. 2006; Fernandez Lorenzo 
et al. 2010), which is in good agreement with our model (but 
see Puech et al. (2008) who find 0.6 magnitudes of bright- 
ening between z ~ 0.6 and 2 = 0). For the B-band TF 
relation, a wide range of evolution has been reported, but 
almost all studies find a brightening in B-band luminosities 
at higher redshifts (e.g., Vogt et al. 1996,1997; Simard & 
Pritchet 1998; Ziegler et al. 2002; Bohm et al. 2004; Bam- 
ford et al. 2005,2006; Weiner et al. 2006b; Chiu et al. 2007; 
Fernandez Lorenzo et al. 2010). Most studies are consistent 
with a brightening of 1±0.5 magnitudes from z = to z ~ 1, 
and thus in agreement with our model. 

Finally, we note that no evolution in the J- or A"-band 
TF relation does not imply no evolution in the stellar mass 
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Figure 10. Evolution of the Tully-Fisher relation using luminosi- 
ties measured in rest-frame optical to near-IR pass bands in our 
disk-galaxy semi-analytic model. At fixed circular velocity, lumi- 
nosities increase out to redshift z ~ 1.5, with stronger evolution 
in bluer passbands than the near-IR. For comparison, the evolu- 
tion of the stellar and baryonic Tully-Fisher relations are shown 
with solid and dotted lines, respectively, which are the same as 
those in the left panel of Fig. [7] 

TF relation. In our model the stellar masses at fixed veloc- 
ity decrease towards higher redshifts, but the light-to-mass 
ratios increase monotonically towards higher redshifts (be- 
cause the stellar populations are progressively younger). The 
net effect is very weak evolution in J- and A"-band TF rela- 
tions between z ~ 2 and 2 = 0. 



4 SUMMARY 

We study the evolution in the zero points of the relations 
between maximum rotation velocity, V, stellar mass, M, and 
rest-frame optical disk size, R, of disk-galaxies in the context 
of ACDM based galaxy formation models. 

Using data from the DEEP2 survey to measure the 
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evolution of the RM relation since z = 1.2, together 
with published results from Conselice et al. (2005), Tru- 
jillo et al. (2006), Kassin et al. (2007), and the SINS sur- 
vey (Cresci et al. 2009; Forster-Schreiber et al. 2009), we 
show that there is a consistent observational picture, with 
one exception, for the evolution of the VMR relations from 
redshifts z ~ 2.2 to z ~ 0. The exception is that the Ha ex- 
ponential disk scale lengths of galaxies from the SINS survey 
measured by Cresci et al. (2009), appear to be a factor of 
~ 3 higher at fixed M than found by other observations. 
This apparent discrepancy can be traced to two factors. 

Firstly, Bouche et al. (2007) and Cresci et al. (2009) 
measure HWHM (half-width half-max) sizes and interpret 
these as exponential disk scale lengths. However, the Ha 
half-light radii from Forster-Schreiber et al. (2009) of the 
same galaxies do not support this interpretation. Instead, 
they suggest that the HWHM overestimates the disk scale 
length by a factor of ~ 2. 

Secondly, using a ACDM based disk galaxy formation 
model (Dutton & van den Bosch 2009) we show that disk- 
galaxies at redshifts z = — 3 are expected to have half-mass 
radii of recently formed stars, J?sf, a factor of ~ 2 higher 
than the half- mass radius of the total stellar mass, R s ta,v 
This is a direct consequence of inside out disk growth. In 
terms of optical half-light radii, our models predict that, 
Rsf — lARi ~ 1.257iV with little dependence on redshift. 

Additionally, since the SINS survey is not a volume lim- 
ited sample of galaxies at high-redshift, and at least some 
of the SINS sample was selected on the basis of disky mor- 
phology with spatially resolved velocity gradients, another 
possibility (which we do not invoke here) is that a selection 
bias against small galaxies exists in the SINS survey. In or- 
der to rule-out this possibility the size- mass relation of SINS 
galaxies needs to be compared to that of a volume limited 
control sample, with the same methods used to derive sizes 
and masses. 

We further show that the observed evolution of the 
VMR relations is consistent with a simple ACDM-based 
model of disks growing inside evolving NFW dark-matter 
haloes. This model adopts a constant disk-to-halo mass 
fraction of md = 0.04, and a median spin parameter of 
A = 0.035, independent of redshift. The galaxy mass fraction 
is consistent with observations at low redshifts (Hoekstra 
et al. 2005; Dutton et al. 2010b), and the spin parameter 
is consistent with expectations from cosmological simula- 
tions (e.g., Bullock et al. 2001b; Maccio et al. 2007). While 
our model is certainly an over-simplification of disk-galaxy 
evolution, it demonstrates that there is no need to invoke 
abnormally high spin parameters to explain the scaling rela- 
tions of star-forming disk-galaxies at high redshifts (z ~ 2) 
as has been claimed by Bouche et al. (2007) and Burkert 
et al. (2009). 

The weak evolution of the galaxy scaling relations since 
z ~ 1 in our model is due to the fact that individual galax- 
ies grow roughly along the scaling relations, and not due 
to weak evolution in the properties of individual galaxies 
themselves. Similar conclusions for the stellar mass- velocity 
relation have been reached by Portinari & Sommer-Larsen 
(2007) using cosmological simulations, and for the size- 
stellar mass relation by Firmani & Avila-Reese (2009) us- 
ing a semi-analytic model similar to the one used here. 
For example, for a galaxy with present day stellar mass of 



3 x 1O 1O M , since z = 1, its stellar mass has increased by a 
factor of ~ 2.5, its half-light size has increased by a factor 
of ~ 2.0, and its maximum circular velocity has increased 
by a factor of ~ 1.15. In terms of samples of galaxies, the 
evolution is just a factor of ~ 1.5 increase in stellar mass, 
at fixed Vmax, and a factor of ~ 1.3 increase in disk size at 
fixed stellar mass, since z = 1. Evolution in stellar mass by a 
factor of ~ 2.5 for galaxies with present day stellar masses of 
3 x 1O 1O M0 is also consistent with the evolution of the star 
formation rate - stellar mass relation since z = 1 (Noeske 
et al. 2007b). 

In our models the evolution of the stellar scaling rela- 
tions is stronger than that of the baryonic scaling relations 
(maximum circular velocity, baryonic mass, baryonic half- 
mass size). This is due to a combination of the inside-out 
nature of stellar disk growth in ACDM cosmologies, coupled 
to a decrease in cold gas fractions with cosmic time. 

In our models the evolution of the baryonic scaling re- 
lations is weaker than that of the virial scaling relations 
of dark matter haloes, assuming a constant galaxy mass 
fraction. For example, the baryonic TF and baryonic size- 
velocity relations evolve by just ~ 0.11 dex from z = 1 to 
z = 0, whereas the corresponding relations between halo 
virial quantities evolve by ~ 0.33 dex from z = 1 to z = 0. 
This difference can be understood as a consequence of the 
ratio between maximum circular velocity and virial circu- 
lar velocity, Kn ax /Vvir, increasing towards lower redshifts. 
This in turn is largely a consequence of the increase in halo 
concentrations with time (e.g., Bullock et al. 2001a). 

While we have shown that there is a consistent observa- 
tional and theoretical picture for the evolution of the VMR 
relations out to redshift z ~ 2, there is much room for 
progress. Our theoretical model makes a number of simpli- 
fying assumptions which are unlikely to be correct in detail. 
We assume that the galaxy mass fraction and galaxy spin pa- 
rameters do not evolve with time. While this is a reasonable 
assumption to start with, we note that this is not a natural 
outcome of our model when we include cooling and outflows. 
Our model also assumes that galaxy disks are smooth, and 
100% supported by rotation. This assumption is valid in the 
local universe, and perhaps up to redshift z ~ 1, but may 
break down at higher (z > 2) redshifts, especially for lower 
mass galaxies. 

Current observational samples at z > 1 are small and/or 
subject to measurement uncertainties and selection biases. 
For the RM relation, it will be possible, in the near future, 
to measure robust rest-frame optical sizes at z < 2.5 using 
large NIR surveys with HST/Wide Field Camera 3. Measur- 
ing maximum circular velocities is currently a challenge at 
high redshifts, as it is limited to ground based NIR spec- 
troscopy. These observations are typically seeing limited, 
which complicates the measurement of maximum rotation 
velocities. However, coming generations of adaptive optics- 
assisted ground based telescopes will open up the field in 
concert with the Atacama Large Millimeter/submillimeter 
Array (ALMA) and the Square Kilometer Array (SKA). 
With this suite of frontier instrumentation, it will be pos- 
sible to measure resolved rotation curves and gas density 
profiles in molecular and atomic gas out to high redshifts. 
Coupled to observations of stellar masses and sizes, these 
observations will enable the evolution of the baryonic VMR 
scaling relations over a large fraction of cosmic time to be 
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measured. These scaling relations will provide a complete 
set of observational constraints with which to test models of 
disk-galaxy formation. 
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